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1. Introduction 

The observation of type la supernovae pOj tells us that the universe is accelerating, 
while cosmic microwave background experiments [2] point to a spatially flat universe. 
Within the context of general relativity, the acceleration is explained by a form of dark 
energy which accounts for 70% of the cosmic energy content with exotic equation of 
state P ~ —p. It is even possible that the effective equation of state of this dark energy 
component has P < —p (phantom energy) [3]. As an alternative to postulating this 
strange form of unseen energy, various modiflcations of gravity have been proposed. 
Among them probably the most popular is the so-called f{R) gravity [4], described by 
the action 

S = J d^x ^ + , (1.1) 

where R is the Ricci curvature, f{R) is a non-linear function of R, g is the determinant 
of the metric tensor gab, and SvrG = 1 (G being Newton's constant) in our notations. 
Other notations follow Ref. [5]. f{R) gravity comes in three versions: the so-called metric 
formalism, in which the action (11. ip is varied with respect to g°'^ and the field equations 
are of fourth order; the Palatini formalism [6], in which the action (jl.ip is varied with 
respect to both the metric and the connection, yielding second order equations, and in 
which the matter part of the action 5'('"'^**e»') is independent of the connection; and the 
so-called metric-affine gravity version in which also 5'(™""e'') depends explicitly on it [7]. 
The physical motivation for all these versions of modified gravity lies in the desire to 
explain the cosmic acceleration without dark energy (and, especially, to get away from 
the embarassing notion of phantom energy). 

Here we restrict the analysis to the metric formalism, in which the field equations 

are 

f\R)Rab - ^ gab = VaVbf - gab^f + T^b (1-2) 

(a prime denotes differentiation with respect to R). 

While several models of metric f{R) cosmology are found in the literature, most 
works consider specific forms of f{R) which are not motivated by fundamental reasons, 
and are usually chosen on the basis of simplicity or ease of calculation. It is more 
interesting not to make assumptions on the form of f{R) and try instead to understand 
this class of theories in as general a way as possible, without choosing the form of f{R). 
In a previous paper we have shown that theories with /" < are not viable due to 
instabilities in the Ricci scalar [8l[9]. Problems with the weak field limit have also been 
pointed out for many choices of f{R) [lO]. Other issues to be studied include the correct 
sequence of cosmological eras (inflation, radiation era, matter era, present accelerated 
era) and whether smooth transitions between them are possible [11]; the presence of 
ghosts and instabilities [12] ; and the well-posedness of the Cauchy problem [8]. 

Here we provide a description of the phase space for general metric f{R) cosmology, 
i.e., without specifying the form of the function f{R). We focus on the geometry of the 
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phase space, the existence of equilibrium points and their stability, and the phase space 
picture of the equivalent scalar-tensor theory. 

Further motivation for our study comes from scenarios of inflation in the early 
universe employing quadratic corrections to the Einstein-Hilbert Lagrangian, i.e., 
f{R) = R + aR^ [I3]. (The phase space of such theories has already received attention 
in the literature |14|.) Quadratic corrections are motivated by attempts to renormalize 
Einstein's theory [15], and are included as special cases in the following discussion. 

It is well known that f{R) gravity can be mapped into the Einstein conformal 
frame in which the theory is equivalent to a scalar field minimally coupled to the 
Ricci curvature ([IE]; see also [EllIH]). Although this version is completely equivalent 
to "Jordan frame" f{R) gravity, at least at the classical level, there have been much 
confusion and misunderstanding about the mapping into the Einstein frame (see, e.g., the 
recent discussion of [19j) and we prefer to proceed directly without using this conformal 
mapping. 

The plan of this paper is as follows. In Sec. El the geometry of the phase space 
of spatially flat Friedmann-Lemaitre-Robertson- Walker (FLRW) cosmology is studied. 
In Sec. [31 the phase space view is compared with that of the equivalent scalar-tensor 
theory, and effective Lagrangians and Hamiltonians are presented for both cases. Sec. [4] 
contains a discussion and the conclusions. 



2. The phase space of f{R) gravity 



In this section, we study the geometry of the phase space of f{R) gravity. Since we 
study regimes dominated by corrections to Einstein gravity (for example, the late-time 
acceleration of the universe, or the early epoch of inflation in Starobinsky-like scenarios 
[T3]). we omit the matter part of the action (11.11) . Motivated by the recent cosmological 
observations, we restrict the analysis to the spatially flat FLRW line element 



ds^ = -dt^ + o?{t){dx^ + dy^ + dz^) 
in comoving coordinates (t,x,y,z). The field equations then reduce to 

'Rf - f 



1 
3? 



2H + 3H^ = - — 



1 



3HRf" 



f"\Rf + 2HRf" + Rf" + -(/ - Rf) 



(2.1) 



(2.2) 



(2.3) 



where an overdot denotes differentiation with respect to t. We assume that /' > to 
have a positive effective gravitational coupling; furthermore, modified gravity theories 
suffer from violent instabilities unless /" ^ ([20l EH [H [9] — see also [22]) and it seems 
appropriate to exclude the case /" < on a physical basis. The field equations are of 
fourth order in a{t). However, when the curvature index = 0, a only appears in the 
combination H = d/a. Since the Hubble parameter is a cosmological observable, it is 
convenient to adopt it as the dynamical variable; then, the field equations (12. 2p and ( 12.31) 
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are of third order in H. The elimination of a is not possible when the curvature index 
7^ 0, or when a fluid with density p = p{a) is included in the picture. We can describe 
the dynamics in the three-dimensional phase space (^H, R, , then the Hamiltonian 
constraint (12 ■2p implies that the orbits of the solutions lie on an energy surface S and 



there is an expression for R for any given value of the other two variables {H, R). This 
is given explicitly by eq. (12.21) as 

where we assume that /" > 0. It is therefore possible to eliminate the variable R given 
a pair {H, R) . Note that there is only one value of R for each given value of {H, R) . 
This situation is different from a general scalar-tensor theory in which one obtains two 
values (j)±{H, 0) for any value of the pair {H, 0) (0 being the scalar field of gravitational 
origin) because the Hamiltonian constraint is a quadratic equation in in which the 
term 0^ is multiplied by the Brans-Dicke "parameter" uj{(f)) (which in general scalar- 
tensor gravities becomes a function of instead of the constant parameter appearing in 
Brans-Dicke theory [23l [24|). The scalar-tensor equivalent of f{R) gravity has u = 
and no quadratic terms in — see the next section. This situation refiects the fact that 
R {H, R) is single- valued. 

By eliminating the variable R, the orbits of the solutions of eqs. (12.21) and (12.31) 
are confined to the two-dimensional, curved, surface S described by eq. ( 12. 4p in the 
three-dimensional space (^H,R,Rj. By contrast, in general scalar-tensor theories with 

c<j 7^ 0, = (f)±{H,(f)) is a double-valued function and the corresponding surface T,± 
is two-sheeted [24]. As a consequence of the two-dimensional nature of S, there can 
be no chaos in the dynamics. Although the standard Poincare-Bendixson theorem [25] 
only applies to a flat, compact and simply connected two-dimensional phase space, it is 
rather straightforward to prove the absence of chaos in the two-dimensional, two-sheeted 
phase space of scalar-tensor gravity [26], and the proof is easily extended to the simpler 
one-sheeted phase space of f{R) gravity. 

As a consequence of the fact that the phase space is constituted of a single sheet, the 
projections of the orbits onto the {H, R) plane can not intersect one another, contrary 
to what happens in the general scalar-tensor case [24j. In general, there is a forbidden 
region in the phase space, corresponding to < and here called = JF^ U JF2, where 

J-i = [{H,R) : Rf'{R) - f{R) - 6HR{H,R) f"{R) < o} (2.5) 

and 

= {{H,R) : f{R)J'{R), or f"{R) are not defined} . (2.6) 

In addition, we assume that /' > to ensure a positive gravitational coupling and that 
/" > to avoid instabilities, which restricts further the region dynamically allowed to 
the orbits of the solutions. The boundary of the connected subset J-'i of the forbidden 
region is 

B^{{H,R,R):f'{R)H' = 0}, (2.7) 



The phase space view of f{R) gravity 



5 



which defines a curve in the surface S, along which the effective energy density Pe/f 
given by eq. p. 141) below vanishes. 

A difference with respect to general scalar-tensor cosmology stands out: in these 
theories, the quadratic equation for the time derivative of the Brans-Dicke-like scalar 
(f){t) provides two distinct roots (p± which differ by a square root term entering with 
positive or negative sign, and there is a forbidden region corresponding to the argument 
of this square root being negative [24j. As a result, the forbidden region intersects the 
two-sheeted energy surface S, creating some "holes" in this surface. The orbits of the 
solutions can change sheet at the boundary of these holes. In f{R) cosmology instead, 
there is no such square root term and the forbidden region does not intersect the energy 
surface S. 

Having chosen {H, R) as variables, the equilibrium points are given by {H, R) = 
(0, 0) and are necessarily the de Sitter spaces 



(iJ, R, r) = {Ho, 12Hl 0) = |^±y A , ^ , j , (2.8) 

where /o = f{Ro), /o = f'{Ro), etc. These fixed points only exist when the condition 
/o//o ^ is satisfied. There is only one condition for the existence of the de Sitter space, 
contrary to general scalar tensor gravity [2lj because eq. (12.21) is not independent of 
eq. (12. 3p and there is only one independent equation left. According to the observational 
data available [I], the present state of the universe is characterized by an effective 
equation of state P ~ —p, i.e., if the f{R) model applies, we are near a de Sitter 
fixed point. 

If /o ^ and /q > 0, two de Sitter equilibrium points (one contracting and one 
expanding) always exist in the R = plane with Hq"^^ = ±./-M-; they degenerate into a 



single point, a Minkowski space, if f{Ro) = 0, in which case they lie on the boundary B 
of the forbidden region J-'i. The latter approaches the energy surface R = R{H,R) at 
these Minkowski space points, which are the only possible global Minkowski solutions. 

The stability of the de Sitter equilibrium points of modified gravity with respect to 
inhomogeneous (space and time-dependent) perturbations was studied in Refs. [27] by 
using the covariant and gauge-invariant formalism of Bardeen-Ellis-Bruni-Hwang |28], in 
Hwang's version valid for generalized gravity [29]. The result is the covariant and gauge- 
invariant condition for linear stability of de Sitter space with respect to inhomogeneous 
perturbations 

^ . (2.9) 

Jo/o 

The contracting de Sitter spaces with Hq < are always unstable [27]. To linear 
order, expanding de Sitter spaces are always stable with respect to tensor mode^. The 
equilibrium points (jT, R, Rj = {Hq, 12Hq, 0) always lie in the {H, R) plane of the phase 
space. 

I It turns out that l|2.9p coincides with the stabiHty condition of de Sitter space with respect to 
homogeneous perturbations [30l I31| . 
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Examples of the phase space geometry for two choices of the function f{R) in the 
gravitational action (ll.ip are presented in fig. [T] and fig. [2l 




Figure 1. Phase space portrait for the model f{R) = R + aR^, with a = 0.0001. The 
vertical axis shows R and R, H, R, a are measured in arbitrary units. 

It is useful to obtain the effective Lagrangian L and Hamiltonian E for f{R) 
cosmology in FLRW space with a and R as Lagrangian coordinates. To the best of 
our knowledge, the effective Lagrangian and Hamiltonian (12.101) and (12. lip below have 
not been presented in the literature, and will be used in future work. For example, they 
can be used in the search for point-like symmetries (see Ref. [32]). An analysis of the 
field equations leads one to conclude that the effective Lagrangian is 

L (a, R, a, i?) = [QH^f + 6Hf"R + f'R-f (2.10) 
and that the effective energ}j§| is 

E (^a, R, a, i?) = (QH^f + QHf'R - f'R + /) . (2.11) 
The canonical momenta conjugated to the generalized coordinates a and R are 

Pa = ^ = 6a^ {2Hf' + f'R) , Pn=^ = Qci'Hf" . (2.12) 

The Hamiltonian constraint ( 12. 2p corresponds to E = 0. By choosing a and R as 
canonical variables, the Euler-Lagrange equation ^ \q^) ~ ^R ~ ^ gives the well 

§ The quantity E can also be regarded as an effective Hamiltonian but it does not contain explicitly 
the generalized momenta. 
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Figure 2. Phase space for the model f{R) — R~ with /x^ = 0.0001, in arbitrary 
units, with R on the vertical axis. 

known relation R = 6 (^^ + while the second equation ^ (^) ^ ^ ~ ^ yields 
eq. (12.31) ■ For non-spatially flat FLRW universes (A; 7^ 0), the Hamiltonian constraint 
takes the form 

k 



— ^ - mkf" 



(2.13) 

a 

in which the scale factor a appears explicitly, and one can not use H as variable instead 
of a. The surface S studied for = separates the region of the phase space accessible 
to the orbits of the solutions for > from the region corresponding to < 0. The 
orbits can not cross the surface S and, in the case = 0, instead, they lie completely 
in S. A similar situation was recognized in an early study of inflation with a massive 
scalar held in the context of Einstein gravity [33] . 



2.1. The effective equation of state 

Using eqs. f l2.2p and (12. 3p and identifying the effective gravitational coupling with 
l/f'{R), one can write the effective energy density and pressure for f{R) models as 
follows: 

Peff = ^^^-3HRf" , (2.14) 
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(2.15) 



The effective energy density pe// is non-negative, as can be seen from inspection of 
eq. (12. 2p . The effective equation of state parameter Weff can be expressed as 



^eff 



P. 



eff 



R^r + 2HRr + Rf" + -(/ - Rf) 



Peff 



Rf'-f 



(2.16) 



- 3HRf" 



Equation (12. 2p guarantees that the denominator on the right hand side of eq. (I2.16P is 
strictly positive, hence, the sign of the effective equation of state is determined by the 
numerator. Type la supernovae yield an effective equation of state w ~ — 1 at present. 
For the model to mimic the de Sitter equation of state with Wejf = —1, it must be 

_ RH -R 

It is convenient to introduce the quantity il){R) = f'{R) to write 



(2.17) 



eff 



-1 + 2 



-1 + 



Rip- f - QHtp 3^/f2 
Alternatively, the deviation from the de Sitter equation of state w 
parametrized by 



Peff + Peff 



ip d 
ip dt 




(2.18) 
-1 can be 

(2.19) 



According to eq. (12.19p . an exact de Sitter solution corresponds to ip = f"{R)R = 0, 
or to ^Z' = Ca{t) = Cooe^"*, where C 7^ is an integration constant. It is easy to 
see that the second solution for ip{t) is not acceptable because it leads to the absurd 
equation f"{R)R = Caoe^"* in which the left hand side is time-independent (for a de 
Sitter solution) while the right hand side is not. 



3. The equivalent scalar-tensor theory 

A modified gravitational action of the form (II. ip can be recast into the form of an 
equivalent scalar-tensor theory ( |34l [35] ; see |16| for a mapping into general relativity 
with a minimally coupled scalar). When /" 7^ 0, the modified gravity action ( ll.ip can 
be rewritten as the scalar-tensor theory 

S=\Jd'xV^ - Vm + , (3.1) 

where 

^(0) = /'(0) , V{<P) = <Pf'-f, (3.2) 

and the Brans-Dicke parameter is = [HH [35]. The Ricci curvature scalar R is 
identified with the scalar field degree of freedom and the mapping to scalar-tensor 
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variables can be seen as a Legendre transformation (see Ref. [I7] for a discussion, and 
Ref. [36] about the physical equivalence). In fact, variation of the action (13.11) with 
respect to yields 

dih dV 

which in turn yields = i? if /" 7^ 0. There is no kinetic term for the scalar in the 
action (13.11) . but this quantity is dynamical because R = (j) obeys the dynamical equation 

In the metric (12. ip the field equations become 

1 r.... ^3_^^ 



V{(f)) - QH^\ 



if = J- ( H^'(f) - ^'(j) - ^"02 I ^ (3 5) 



and ( 13. 3p . or 



1 

24) 



H=^, [Hf"4> - /"0 - /'"02 j ^ (3.7) 

{R -<l>)f = 0. (3.8) 

By viewing the Hamiltonian constraint (13. 4p as an algebraic equation for 0, one can 
eliminate from the three-dimensional phase space (^H, 0, 0^ (which corresponds to the 
phase space (if, R, R) of Sec. [2|), obtaining 

{H, 0) = ^ (0/ - / - e/'ii^) . (3.9) 

This represents a one-sheeted energy surface, which corresponds to the special case u = 
of scalar-tensor cosmology |24| . as remarked in the previous section. The equilibrium 
points (ii, 0) = (0,0) are de Sitter spaces with constant scalar field and satisfy 



(^o,0o)= |±W^,^j ; (3-10) 

they lie in the = plane. In general scalar-tensor gravity it is possible to have de 
Sitter solutions with non-constant scalar field (which, however, are not fixed points). In 
the scalar-tensor equivalent of f{R) gravity this is not possible because of the constraint 
(p = R = 12Hq which forces to be constant for all de Sitter solutions, which have 
constant Hubble parameter Hq. 

The gauge invariant stability condition for these de Sitter equilibrium points with 
respect to inhomogeneous perturbations, obtained in Ref. [30], is again 



ifoLjplK^O^ (3.11) 
Jo Jo 
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The effective Lagrangian and Hamiltonian for the equivalent scalar-tensor theory (13 .ip 
are now given by 



These can be obtained from well known Lagrangians of scalar-tensor theories ([37] and 
references therein) and do not constitute an original result, contrary to the L and E 
given by eqs. (l2T0D and fl2TTD . 

4. Discussion and conclusions 

The phase space view of a spatially homogeneous and isotropic universe is useful for 
building models of the current acceleration of our cosmos and in trying to discover what 
fuels it. The high degree of symmetry of the FLRW metric (12.11) considered turns the 
Einstein equations into ordinary differential equations, for which the Cauchy problem 
is trivial and is covered by well-known theorems. Instead, the initial value problem of 
f{R) gravity for general metrics has not been discussed, except for specific choices of 
the function f{R) [Ml El], and will be addressed elsewhere. 

The dynamical system approach is a powerful tool which is appropriate to study the 
dynamics for general initial conditions when exact solutions can not be obtained (i.e., in 
most situations). As such, dynamical systems theory is widely applied in cosmology (see, 
e.g., Refs. [HUSO]). The theory of dynamical systems has been applied repeatedly also to 
f{R) cosmology, but usually for special choices of the function f{R) [41j- Our goal is to 
study the phase space of f{R) cosmology without committing to a specific choice of the 
function f{R). Motivated by the recent observations of type la supernovae [I] and by 
the current cosmic microwave background experiments [2], we have restricted the scope 
of our analysis to spatially flat FLRW cosmologies. Advantages of our approach include 
its generality, the use of physical observables as dynamical variables, and the ease of 
comparison with the general picture of the phase space for the equivalent scalar-tensor 
gravity. The obvious limitation is that the phase space portrait is necessarily incomplete 
without choosing a specific form of /(i?). However, in the absence of compelling physical 
indications on the form of this function, it is hoped that our results on the geometry of 
the phase space, the equilibrium points and their stability, and the effective Lagrangians 
and Hamiltonians can be used as a preliminary step and a guide to further understanding 
the cosmological dynamics in any specific f{R) theory. The extension to more general 
gravity theories of the form / (i?, RabR"^'', Rahcd.R"'^'^'^) will be discussed elsewhere. 
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